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A petroleum reservoir is a porous medium that contains hydrocarbons. The major goal 
of any reservoir study is to predict future performance of the reservoir and find ways and 
means of enhancing the ultimate recovery under various operating conditions. 
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Classical reservoir engineering deals with the reservoir on a gross average basis (tank 
model) and cannot account adequately for the variations in reservoir and fluid parameters 
in space and time. 



Classical Reservoir Engineering Methods: 

• The analogical methods utilize features of mature reservoirs that are analogous to 
the target reservoir in an attempt to forecast the performance of this target zone or 
reservoir. 

• The experimental methods measure physical properties, such as pressure, 
saturation, and/or rates, in laboratory cores and then scale them up to the whole 
hydrocarbon accumulation. 

• The mathematical methods use model equations to forecast reservoir performance. 
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The mathematical methods 

Here we briefly mention the mathematical methods because of their direct relationships 
with reservoir simulation. 
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1. Material Balance Methods: 

The classical material balance methods use a mathematical representation of a reservoir 
or drainage volume. Their basic principle is based on mass conservation; i.e., the amount 
of mass for water, oil, or gas remaining in the reservoir after a production period equals 
the difference of the amount of mass originally in place and that was removed from the 
reservoir due to production, plus the amount of mass added due to injection and 
encroachment. 

2. Decline Curve Methods: 

The classical decline curve methods use one of three mathematical declines (exponential, 
hyperbolic, and harmonic) to describe the rate of oil production decline. 

= -if, 

q at 

Historical data are matched by choosing the parameters C and b that minimize the error 
(often using the least squares error) between the data and the equation. Then 
extrapolating the historical data into the future to predict reservoir performance using the 
matched equation. A major assumption of any extrapolating method is that all processes 
occurring in the past will continue in the future. In addition to not giving any spatial or 
temporal information the decline curve methods cannot be used for “what-if” scenarios. 
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Statistical methods employ empirical correlations that are statistically obtained using the 
past performance of some reservoirs to forecast the future performance of others. They 
are a generalization of the analogical methods. A correlation is developed with data from 
mature reservoirs in the same region, with the same lithology (e.g., carbonate or 
sandstone) and under the same operating conditions (e.g., solution gas drive or water 
flood). For the reservoir engineer to be confident in using an empirical correlation model, 
reservoir properties must be within the limit of the regression database used to develop 
such a model. 



4. Analytical Methods: 

Analytical methods, such as pressure-transient and Buckley-Leverett methods, use the 
analytical solution of a mathematical model. The model consists of a set of differential 
equations that describe the flow and transport of fluids in a petroleum reservoir, together 
with an appropriate set of boundary and/or initial conditions. To solve these equations 
exactly, simplifying assumptions must be made to reduce the complexity of the model. In 
general, these assumptions are very restrictive. For example, in the pressure-transient 
method, the assumptions require single-phase flow in a horizontal reservoir with uniform 
thickness and small pressure differences under the laminar flow condition. The Buckley- 
Leverett method for two-phase flow ignores gravity and capillary forces under the 
incompressibility condition. However, since much of the physics of a problem is 
preserved, the analytical methods are often used to determine how various parameters 
influence reservoir performance. Furthermore, these methods can be used to validate 
reservoir simulators. 
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The physical model derived for identifying a petroleum reservoir is very complex and 
contain many unwanted effects of uncertainties and disturbances since the reservoir 
is dynamic model which continuously change in nonlinear behavior. In recent 
publications there are some attempts to replace large scale reservoir models with 
simplified proxy models such as Capacitance Resistive Models. However there are 
another techniques are implemented lately such as Fully data driven black box models. 
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L. P. Dake, The Practice of Reservoir Engineering, 

“There should be no competition between material balance and simulation, instead they must be 
supportive of one another: the former defining the system which is then used as input to the 
model. Material balance is excellent at history matching production performance but has 
considerable disadvantages when it comes to prediction, which is the domain of numerical 
simulation modelling. ” 



Material balance has been used in the industry for the following main purposes: 

1. Determining the initial hydrocarbon in place (e.g. STOMP) by analyzing mean reservoir 
pressure vs. production data; 

2. Calculating water influx i.e. the degree to which a natural aquifer is supporting the production 
(and hence slowing down the pressure decline); 

3. Predicting mean reservoir pressure in the future, if a good match of the early pressure decline 
is achieved and the correct reservoir recovery mechanism has been identified. 

4. Predicting the reservoir performance by introducing time-flow rate relationship. 
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Formation Volume Factor in the Black Oil model 



The formation volume factors (FVF) tell how much the 
oil, gas and water is compressed at a given pressure. 

B 0 = reservoir volume of oil / standard volume of oil 

B g = reservoir volume of gas / standard volume of gas 

B w = reservoir volume of water / standard volume of 
water 



B 0 vs. P 




Bg vs. P 




B w vs. P 
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Solution Gas-Oil Ratio in the Black Oil model 



The Rso plot shows how the solution gas ratio develops 
vs pressure. When the pressure reaches the bubble 
point pressure, it is no longer possible to solve more gas 
into the oil. Thus the gradient of the curve becomes 
zero. 

R s = standard volume gas / standard volume oil 



Rso vs. p 




P 
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Material Balance Methods: 



Symbols 


Units 


Units SI 


B, 


Gas formation volume factor 


bbl/SCF 


MVSCM 


Ef 


Oil formation volume factor 


bbl/STB 


M 3 /SCM 




Total formation volume factor 


bbl/STB 


M 3 /SCM 


B w 


Water formation volume factor 


bbl/STB 


M 3 /SCM 


c f 


Pore compressibility 


vol/vol/psi 


vol/vol/Mpa 


c ™ 


Water compressibility 


vol/vol/psi 


vol/vol/Mpa 


G 


Initial gas -cap volume 


SCF 


SCM 


G 

P 


Cumulative gas produced = G pi + 


SCF 


SCM 


% 


Cumulative solution gas produced 


SCF 


SCM 


% 


Cumulative gas cap produced 


SCF 


SCM 


m 


Ratio initial reservoir free gas volume to 
initial reservoir oil volume 


bbl/bbl 


M 3 /M 3 
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N 


Stock tank oil initially in place 


STB 


STM 3 


N 

P 


Cumulative tank oil produced 


STB 


STM 3 


P 


Average reservoir pressure 


psi 


MPa 


Pl 


Initial reservoir pressure 


psi 


MPa 


R 

■n 


Cumulative gas/oil ratio 


SCF/STB 


SCM/STM 3 


t 

R 

£, 


Solution gas/oil ratio 


SCF/STB 


SCM/STM 3 


S 

TV 


Average connate water saturation 


fraction 


fraction 


W 

e 


Cumulative w'ater influx 


bbl or STB 


M 3 or STM 3 


w p 


Cumulative w'ater production 


bbl or STB 


M 3 or STM 3 




Other subscripts 






i 


at initial conditions 






b 


at bubble point 
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Material Balance Methods: 



Amount of fluids present 




Amount of 




Amount of fluids remining 


in the reservoir initially 


> — < 


fluids produced 


> = < 


in the reservoir finally > 


(st. vol.) 




(st. vol.) 




(st. vol.) 



Note that “fluids produced” include all influence on the reservoir: 

• Production 

• Injection 

• Aquifer influx 



on 

producer 



Gas Water 

Injector Injector 
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Material Balance Methods: 



MB for Oil Reservoirs Above Bubble Point: 
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Material Balance Methods: 



MB for Oil Reservoirs Below Bubble Point: 
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MB with gas cap and water drive: 




Pi 



P 




(W e - W p ) + (G-G pH; )B g -GB^ 

= NB oi -((N-N p )B 0 +(NR sj -(N-N p )R s -G ps )b„) 
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Effect of Pore Volume Changes: 

1 . Impact of pore volume changes due to rock: 



Overburden pressure 
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Material Balance Methods: 



Effect of Pore Volume Changes: 

1 . Impact of pore volume changes due to connate water: 



Compressibility of water: 



1 

p 


AV pw 


w V 

pw 


Ap 



AV = c ApS V 

pw w r wc p 



Total Pore Volume Change due to rock & water: 



AV = AV +AV 

p pr pw 

.-.AV =(c f +c S )ApV 

p V f w wc ) r p 
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Net water influx + gas cap expansion + pore volume reduction = 



/ 1 



\ 



0 W. - W p B w ) + (G - G„ ) B s - GB,, + A pNB cl 




Original oil volume - volume of remaining oil and free solution gas. 
I 
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v Material balance in the form of straight line : 




Material balance not a difficult concept but it is difficult in applying it to real reservoirs. There is 
often inadequate understanding of drive mechanisms. 




Odeh & Havlena (1963) rearranged MB equation into different linear forms. Their method requires 
the plotting of a variable group against another variable group selected depending on the drive 
mechanism. 




The technique is referred to as history matching. 


| faq| [references about 


HELP 




Reservoir Simulation 



• INTRODUCTION 

Classical Methods 
Mathematical Methods 
Material Balance 

MODELLING 

APPLICATIONS 

TUTORIALS 



Material Balance Methods: 

Predicting the reservoir performance using Material Balance Equation: 
Material balance in the form of straight line : 

Material balance can be rewritten as: 



- w i.j 

= N[(B a -Bj+(R„-R ! )B s \ 
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+ mNB oi 



B 
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V s , 
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Material Balance Methods: 



Predicting the reservoir performance using Material Balance Equation: 
Material balance in the form of straight line : 



Havlena and Odeh simplified equation to:- 



F = NE 0 +NmE g +NE fw +W e 

F = N p [B„ + ( R p - R s )B t \ 

E„ =(B„ -Bj + iR, -fiJS s ....bbl/STB 



Efw ~ 
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Predicting the reservoir performance using Material Balance Equation: 
S Material balance in the form of straight line : 



Tarek Ahmed, Reservoir Engineering Handbook, Third edition, Chapter 11 

1 . Case 1 : Determination of N in volumetric undersaturated reservoirs 

2. Case 2: Determination of N in volumetric saturated reservoirs 

3. Case 3: Determination of N and m in gas cap drive reservoirs 

4. Case 4: Determination of N and We” in water drive reservoirs 

5. Case 5: Determination of N, m, and We in combination drive reservoirs 

6. Case 6: Determination of average reservoir pressure, p 
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Material Balance Methods: 



The particular advantage of material 
balance models is that they are very 
simple. They can address questions 
relating to average field pressure for 
given quantities of oil/water/gas 
production and water influx from given 
initial quantities and initial pressure 
(within certain assumptions). 



(a) 'Tank 1 ' Model of the Reservoir 



Wells Offtake 



Average Pressure 
Average Saturations 



P 

S o , =£ w and S 




Aquifer 
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However, the material balance model is 
essentially a tank model, it cannot 
address questions about why the 
pressures in two sectors of the reservoir 
are different (since a single average 
pressure in the system is a core 
assumption). 



The sector model is somewhat more 
complex in that it recognizes different 
regions of the reservoir. 



(b) Simple Sector Model 





Producer - West Flank 


Producer - East Flank 








— 
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Oil Leg 

4 




Aquifer 
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Material balance can not answer the 
more detailed questions especially for 
mature field like: 

where should I locate an infill well and 
should it be vertical, slanted or horizontal 
? 



For such complicated questions, the 
model has to be more detailed and it 
contains more spatial information. Like 

Simulation Model 



(c) Fine Grid Simulation Model of a Waterflood 
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A simulation model shows the main features of 
a real system or resembles it in its behavior , but 
is simple enough to make calculations on. These 
calculations may be "analytical" or "numerical". 
By analytical it means the equations that 
represent the model can be solved using 
mathematical techniques such as those used to 
solve algebraic or differential equations. An 
analytic solution would normally be written in 
terms of “well known” equations or functions (x2, 
sin x, ex etc).. 



Learning goals 

• Basic understanding of material balance 
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